Efficiency of structured adiabatic quantum computation 
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We show enough evidence that a structured version of Adiabatic Quantum Computation (AQC) 
is efficient for most satisfiability problems. More precisely, when the success probability is fixed 
beforehand, the computational resources grow subexponentially in the number of qubits. Our study 
focuses on random satisfiability and exact cover problems, developing a multi-step algorithm that 
solves clauses one by one. Relating the computational cost to classical properties of the problem, 
we collect significant statistics with up to N — 140 qubits, around the phase transitions, which is 
where the hardest problems appear. 



While the general framework of Quantum Computa- 
tion is to a large extent justified by the progressive minia- 
turization of integrated circuits, there is still a major 
open problem in the field which is the finding of new 
quantum algorithms. This task is rather difficult due 
to both the lack of a fully quantum programming model 
and the general requirement that any quantum algorithm 
should outperform its classical counterparts -a perhaps 
too ambitious expectation which imposes both proofs of 
the classical and quantum complexities. 

The Adiabatic Quantum Computation paradigm [l[ is 
a very appealing framework for the development of new 
algorithms. While equivalent to the circuit model 
this paradigm aims at adiabatically preparing the ground 
state of a quantum Hamiltonian, a state which encodes 
the solution to a problem we want to solve. The re- 
source that measures the computational cost in AQC is 
the time required to prepare the state and the adiabatic 
theorem Q provides a criterion to bound this time, ei- 
ther numerically or analytically. In addition to this, the 
adiabatic formulation is particularly well suited for con- 
straint solving problems. The constraints may be en- 
coded as energy penalties in a Hamiltonian, so that any 
solution of the whole problem must be a ground state of 
the final Hamiltonian. 

The traditional interest on AQC is precisely rooted on 
its relation to constraint solving and in particular to logi- 
cal satisfiability (SAT) problems. From a practical point 
of view, many real- world problems, including automated 
hardward design and verification, can be suitably repre- 
sented in SAT form. Furthermore, the subset of 3SAT 
problems were also the first problems shown to be NP- 
complete [J, [5| . This implies that if we had an algorithm 
to efficiently solve 3SAT, we could also solve all problems 
in the much larger NP family [28(. However, to date, 
the extensive literature about classical SAT algorithms 
suggests that these problems are exponentially hard to 
solve @] , even in the average and typical cases 0| ■ 

Regarding the utility of quantum computation in solv- 
ing SAT problems, there is not yet a clear cut answer. De- 
spite the promising results of the very first AQC simula- 
tions [H, nowadays AQC is expected to be expontentially 



costly for the worst-case instances of NP-complete SAT 
problems no use of the structure of 

the problem is made. One may argue that the worst- 
case complexity is less relevant than the average-case 
cost [HI], and there are evidences for average exponen- 
tial speedups in the quantum query problem [l3| . Many 
works have adopted this point of view and studied the 
running times of avera ge or t y pical adiabatic compu- 
tations 0, [3, [H, [H, [Hi . Here the results are 
still mixed, ranging from polynomial [l8| to subexpo- 
nential [H, [l6| and to exponential [3] growth. These 
discrepancies can be attributed to the difficulty of simu- 
lating a quantum algorithm, which limits the analysis to 
small instances and statistical samples. 

In this work we want to rule out not only the worst-case 
restriction, but also the use of unstructured algorithms. 
We will introduce a new composite algorithm for SAT 
that alternates diabatic and adiabatic stages. We will 
then relate the running times of the algorithm to classical 
properties of the SAT instances. Using this connection, 
we exactly characterize the performance of the algorithm 
for very large problems, with a much larger sampling 
space than any previous study — up to N = 140 bits 
with 2 x 10 6 instances — . With this data we will answer 
a question of great practical relevance: given a desired 
success probability, what is the scaling of the resources 
required to solve that fraction of randomly picked, but 
classically hard problems. The result will be that the re- 
sources grow subexponentially, even in regimes in which 
the classical algorithms would be exponentially costly 0] • 

This work focuses on the solution of SAT problems. 
An instance of this family is defined by a set of M logical 
constraints or clauses on a finite set of N Boolean vari- 
ables (si, S2 ■ ■ ■ sn). In the particular case of p-SAT, each 
clause is a Boolean function depending only on p variables 
f k (s kl ,s kz ,...s kp ) : {0,1}®p^ {0,1}. The problem is to 
determine whether there exists an assignment of the N 
variables satisfying all M clauses, a task for which there 
exist both heuristic and probabilistic classical algorithms. 

Our quantum algorithm for SAT is a hybrid that alter- 
nates non-adiabatic evolution with adiabatic steps, in a 
way that goes beyond structured AQC [T3|. The key in- 
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FIG. 1: Efficiency of the quantum algorithm. In (a) we plot 
the total running time vs. the number of qubits. We have used 
Eq. |T| with n = 1. We plot average results for X3SAT (cir- 
cles) and 3SAT (crosses) together with standard deviations. 
The inset shows the deviation from a linear law for X3SAT. 
Figure (b) is a log-log plot of the average value (solid) and 
variance (dashed) of the value /3, which is the logarithm of 
the maximum time required by one of the adiabatic intervals. 
The lines are fits to power laws, /3, A/3 ~ aN b . Each point in 
the plots involves 2 x 10 6 random instances. 

gredients are to initially sort the set of logical clauses, and 
to have each adiabatic step solving a single clause. Let 
us call S m the set of assignments which are compatible 
with the m first logical clauses. As shown below, during 
the m-th time interval of adiabatic evolution (t m _i,i m ), 
the quantum register will evolve from a linear superpo- 
sition of all states in S m -\ to a linear superposition of 
those in S m . For this step of the computation to succeed 
we require a run time 

tm ~ trn — l = £^ X d m ^\ j d m , (1) 

where d rn = \ S m | is the size of the set S m and fl is some 
big multiplicative constant specified by the adiabatic the- 
orem [3J. This relates the resource of the adiabatic algo- 



rithm — time — to a classical property — the structure of 
solution spaces — , and it will be the main tool in our anal- 
ysis of the efficiency. Note that if we add all clauses in 
the same adiabatic interval, the adiabatic time becomes 
exponentially large T soarc h ~ fMo /c?m ~ 2 , consiste ntly 
with search algorithms for unstructured problems [201 ] . 
Let us now present the algorithm and derive Eq. (TT]). 

The algorithm. In AQC, the solution of a compu- 
tational problem is encoded in a state ip g which is the 
ground state of a problem Hamiltonian, Hp. To prepare 
tjjg, the Hamiltonian of the quantum register is slowly 
modulated from an initial operator H(0) with an eas- 
ily prepared ground state ip{0), up to our final operator 
H(T) = Hp. The adiabatic theorem states that if the 
evolution is slow enough, the quantum register will end 
up in the desired state, ip(T) = ip g . The condition "slow 
enough" evolution means that the changes in the Hamil- 
tonian happen in a long time scale T ^> l/(Ai?) 2 , larger 
than the minimal gap AE between the ground state man- 
ifold and the excited states of Hit). 

In our algorithm the initial state is a linear combi- 
nation of all states in the computational basis |^(0)) oc 

EJ S ) = (1°) + \ 1 ))® N ■ This product state is evolved 
with a Hamiltonian that is piecewise continuous. Within 
the m-th interval, t £ (t m _i, t m ) , we write 

H(t) = H^ + [1 - \ m {t)]H { B ] + Kn(t)H^\ (2) 

with adiabatic parameter A m (t) = (t— t m -i)/(t m — t m _i). 
The three operators are 

Ftp 1^) = \}- fm{Smn ^m2i ^3)] |s) , (3) 

H B l) \s) = -^El 8 ')' and W 

s' 

rn— 1 

H ( c m) = r^ifW. ( 5) 

k=l 

Respectively, a term that penalizes invalid configurations 
for the m-th clause, a component that mixes all possible 
configurations and an operator that penalizes any viola- 
tion of the m — 1 clauses we have solved so far. We expect 
that, if the adiabatic steps succeed and the non-adiabatic 
steps do not introduce significant errors, after each time 
t m the register will be in a linear combination of all as- 
signments which are compatible with the m first logical 
clauses, ip{t m ) ~ |S m ) oc J2 s es m I s ) • 

Non-adiabatic errors. In between adiabatic intervals 
the Hamiltonian changes non-adiabatically from ff(t^) 
to iJ(t+ ). This implies strengthening the last clause by a 
factor r and switching on the term H^. Without affect- 
ing the execution time, this process introduces a control- 
lable error on the state, which amounts to the difference 
between the ground states of H(t), that is ipgi^i imme- 
diately before t m and after it. Applying the results in 
Ref. [2| we obtain (^ s (^)|V'fl(*m)) — 1 — r^T- Since the 
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difference is 0(1/F), the M adiabatic steps will amount 
to an overall error M/T, which can be decreased with 
resources that are polynomial in M. 

Cost of the algorithm. Given the previous assump- 
tion that the constraints term is dominant, r 3> 1, 
we may diagonalize the model Hamiltonian |(5J) approxi- 
mately, and thus compute the minimum energy gap dur- 
ing each adiabatic step. We separate the Hilbert space 
into two parts, 7i = TL m -\ © "H-m—ii the s P ace °f so ~ 
lutions which are compatible with the previous to — 1 
clauses, Ti m -i = lin{|s) ,s e S m -i}, and its orthogonal 
complement. During each interval (t m -i, t m ), the low en- 
ergy spectrum of H(t) is approximately given [2J by that 
of its projection onto W m -i, an operator that we call 
H m —i(t). In the two dimensional space spanned by the 
vector S m and its complement |H^) oc Ylses ns I s ) > 
we may write i? m _i = (A„ — 5)!+ §(A — 1) sm(2a n )cr x + 
\ [(An — 1) cos(2a n ) — A] a z , expressed using the Pauli 
matrices and cos(a„) = y/ d n /d n -i = (3 m |H m _i). This 
matrix has a minimal separation or energy gap at A„ — 
1/2, given by AE rn = cos(a TO ) = y/d m /d m -i- This pro- 
vides a first estimate for the running time of the adi- 
abatic step, t m - tm-i > l/(AE m ) 2 or Eq. |T]). But 
for this reasoning to be complete, we must also con- 
sider the influence of states with energy » T. Following 
Ref. 0] the high energy states produce corrections in the 
energy gap and the ground state wavefunctions which 
are small, 0(1/T). Therefore, by imposing a penalty 
r ~ K^/dnZi/dn at each step, with a constant k 3> 1, or 
using a global T larger than all the individual factors, we 
can ensure that Eq. (fT]) still holds. Finally, since multi- 
plicative factors in the Hamiltonian affect the total run- 
ning time, one should actually consider an adimensional 
running time Q such as T = Tmax s ||if(s)|| ~TxT, but 
with the previous choices, both T and T scale similarly. 

Performance. For our benchmarks we have chosen 
two sets of problems: Exact Cover and 3SAT in conjunc- 
tive normal form (CNF). The first family, also known 
as X3SAT (or one-in-three SAT), has clauses that re- 
quire exactly one bit to be set: f k (s kl , s k2 , s ka ) = 1 iff 
Ski + s k2 + s k 3 = F These problems have been exten- 
sively used to test the efficiency of unstructured adia- 
batic quantum computation [U [TEl. lit 21 1. However, it 
is known that even though they form an NP-complete 
set, Exact Cover problems tend to be easier than gen- 
eral 3SAT ones. Hence, we have also studied this more 
general set of problems, casting them in a convenient 
form: f k (s kl ,s k2 ,s ka ) = iff s kl = a k , s k2 = b k , s ka = 
c k , where each function forbids a precise combination 
(a k , b k , c k ) of the bits involved. 

Our statistical study is based on a random sampling 
of 3SAT and X3SAT problems, generated with a bias to- 
wards hard instances. This is done by fixing the ratio of 
clauses to bits, a = M/N, close to the phase transition 
from trivially satisfiable problems to problems with no 



solution. It is known that the density of hard problems 
is highest around these narrow regions [22I I23I [24. 12R |26| . 
as at G (4,4.24) and a X 3SAT S (0.54,0.64), while 
the outer regions can be easily solved using specialized 
solvers. For each instance we generate M clauses, each 
one acting on three different random bits. In the case of 
3SAT we also generate the random bit sequences to be 
rejected, (a k ,b k ,c k ). Our pseudo-random number gen- 
erator is a Mersenne- Twister (2?} which has a period 
~ 2 19937 and good equidistribution properties. 

We have applied our algorithm to a sufficiently large 
selection of randomly generated 3SAT and Exact Cover 
instances. For each random instance of 3SAT or X3SAT 
the clauses are sorted according to the indices of the bits 
that are involved. We loop over the sorted clauses, it- 
eratively creating the space of valid assignments, which 
is decomposed into configurations of the n a (m) active 
and N — n a (m) inactive bits, S m = A m (g> I m . Since 
d m = \A m \2 N ~ na ( m \ only A, m needs to be stored and 
computed. The set A m itself is created from A m _i us- 
ing a parallelized algorithm which involves 128 processors 
and a total of 256Gb of memory for the largest problems. 
Note that sorting the clauses is an essential ingredient to 
make both n a (m) and \A m \ grow slowly. 

The optimal total running time tu is computed as the 
sum over all clauses of the terms in Eq. |T]) . As shown in 
Fig. UK it has an average behavior which is close to linear 
in the number of bits. However, preparing an experiment 
that runs in this optimal time requires some knowledge 
about the solution spaces, as some intervals require more 
time than others. A more realistic goal is choosing a 
running time that ensures a high probability of success for 
arbitrary instances. For this we introduce a new variable 



(3 = maxlog 2 (d m _i/d m ), 



(6) 



with which we can uniformly bound the total running 
time of a given problem as 



t M < QM x 2 







(7) 



We have computed the value (3 for every instance in 
the samples. By studying the statistics of this quan- 
tity over different problems we infer an optimal value, 
i3 op t(o~, N, M), such that running the quantum computer 
for a time T opt = f2M exp(/3 opt ) ensures the correct solu- 
tion of a large fraction a (say a — 99.99%) of all iV-bit 
hard instances. The scaling of (3 op t with respect to the 
problem size and success probability is a more meaning- 
ful characterization of the algorithm than the scaling of 
the worst case times. 

The first results are shown in Fig. [TJa, where we 
plot the average value of (3 and its variance, and fit 
them to curves of the form [3x3SAT = 1.14iV - 162 , 
APxssat = 0.1A - 447 , and (3 3SAT = 0.4687V - 438 , 
AfosAT = 0.1337V - 581 . Based on this, the running time 
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FIG. 2: Probability distribution of algorithm running times. 
The inset shows the probability distributions, P(/3;N), of 
the logarithm of the interval that takes the most time in 
our algorithm @, for N = 8,16,32,60 and 140 bits (left 
to right). These distributions are subject to an exponential 
fit P((3;N) ~ exp(— /?/£), which gives the parameter £(N) 
shown in the plot. This value grows subexponentially in the 
number of qubits and is a realistic measure of the algorithm 
efficiency. 



without previous knowledge of the problem could be 



T opt ~ fi x N x 2 



0+6A0 



N2 



■:N r 



(8) 



where exponent 6A/3 guarantees a very small probability 
of failure. This suggests an overall behavior that follows a 
sub-exponential law with r < 0.5 for X3SAT and r < 0.6 
for 3SAT. 

By studying the actual probability distribution of f3 
we can confirm the sub-exponential behavior, obtain a 
better choice of the running time and actually estimate 
the failure probability. Since (3 cannot be larger than the 
number of bits, these distributions have all finite support. 
Nevertheless, as shown in Fig. [2J their tails can be accu- 
rately bounded by an exponential law exp(— /?/£). Ac- 
cording to this law, the failure probability of a particular 
setup which uses a running time T opt = CIM exp(/3 opt ) is 
bounded by 



Pfail = 



Pop 



(9) 



Our numerical analysis reveals that £ grows algebraically 
with the number of bits. Therefore, a fixed success prob- 
ability can be achieved with subexponential resources. 

Summing up, in this work we have presented a hybrid 
quantum algorithm which has the advantage that its ef- 
ficiency can be computed at a relatively low cost. We 
have shown that it is possible to establish a running time 
such that the algorithm solves a large fraction of ran- 
domly picked but still classically hard problems. This 



time is found to scale subexponentially in the size of the 
problem, which is to be compared with the average ex- 
ponential behavior of classical algorithms 0]. While the 
long range interactions make this algorithm suboptimal 
for implementation, our work opens the path to the de- 
velopment and analysis of other structured algorithms. 
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